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Abstract 

Transverse polarization of 3 Si charmonium states, produced directly in pp 
collisions at asymptotically large transverse momentum p t , has emerged as 
the most prominent test of color octet contributions and spin symmetry in 
quarkonium production. We present predictions for the polar angle distribu- 
tion at moderate values of p t ~ 4 — 20 GeV, covered by the Tevatron Run I 
data. We update the fits of NRQCD matrix elements and discuss their the- 
oretical uncertainties. With our best fit values, no transverse polarization is 
expected at p t ~ 5 GeV, but the angular distribution is predicted to change 
dramatically as p t increases to 20 GeV. 



PACS numbers: 14.40. Gx, 13.88.+e, 13.85.Ni 



During the past three years, the phenomenology of quarkonium produc- 
tion has undergone a phase of exciting developments. These followed from 
the application of Non-Relativistic QCD (NRQCD) ||, an effective theory 
that disentangles physics on the scale of the heavy quark mass tuq, relevant to 
the production of a heavy quark pair, from physics on the scale of the bound 
state's binding energy ttiqv 2 , relevant to the formation of the quarkonium. 
Inclusive quarkonium production is now viewed as a two-step process, where 
the production of a heavy quark pair QQ[n] in a certain angular momentum 
and color state n is described by a perturbatively calculable short-distance 
cross section da n and the subsequent quarkonium formation is parametrized 
by a nonperturbative matrix element (0%), which is subject to the power 
counting rules of NRQCD. 

For J I ip and ip' (collectively denoted by ip) only the color singlet n = 
intermediate state contributes at lowest order in the non-relativistic veloc- 
ity expansion in v 2 . (We use spectroscopic notation with the superscript 
in brackets denoting the color state.) Three more nonperturbative param- 
eters, related to color octet configurations n = 3 s[ 8 \ ^S^, 3 Pj 8 ^ appear as 
v 4 -corrections. Their importance for ip production has been understood in 
[§]. Because gluons couple easily to 3 S\ octet states, gluon fragmentation into 
color octet quark pairs appears as the most plausible explanation of the large 
direct ip production cross section observed at the Tevatron ||. Subsequent 
to this observation many charmonium production processes have been recon- 
sidered to gauge the impact of color octet contributions and to determine the 
corresponding nonperturbative matrix elements. (For reviews and references 
see [|J.) Among the specific predictions of NRQCD, transverse polarization 
of direct ip (i.e. not from B and Xc decays) at large transverse momentum 
Pt in pp collisions at the Tevatron has emerged as the most distinct and 
most accessible signature |(J. Transverse polarization also discriminates 
NRQCD from other approaches, like the color evaporation model which, in 
its most naive version, predicts the quarkonium to be produced unpolarized. 
At moderate p t the polarization signature provides an important consistency 
check on the relative weight of the intermediate color octet states inferred 
from the unpolarized differential cross section. 

Transverse polarization at large pt is a direct consequence of gluon frag- 
mentation as well as spin symmetry of the leading order NRQCD effective 
lagrangian. When p t ^> 2m c the fragmenting gluon is effectively on-shell 
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and transverse. The cc pair in the color octet 3 Si state inherits the gluon's 
transverse polarization and so does the if), because the emission of soft gluons 
during hadronization does not flip the spin. Both, corrections due to spin- 
symmetry breaking and higher order fragmentation contributions have been 
estimated to be rather small ||. However, the color octet ^Sq 8 '* and 3 F'j 8 ' ) 
production channels, which do not have a fragmentation interpretation at 
leading order in the strong coupling a s , have large short-distance coefficients 
[0] and can dominate the pt-distribution over a significant range of pt covered 
by the present Tevatron data. The polarization yield from these subprocesses 
has not so far been calculated. 

In this Letter we report on the calculation of all subprocesses for polarized 
direct if) (not from B and \c decays) production at leading order in a s . We 
update the extraction of NRQCD matrix elements from unpolarized cross 
sections first performed in J7[ and discuss their uncertainties. The prediction 
for the polar angle distribution includes the higher order fragmentation cor- 
rections from |H and covers the entire pt-range pt ~ 4 — 20 GeV accessible 
with the Tevatron run I data. The polarization analysis of Tevatron data is 
under way ||. 

The differential cross section for pp — > ip^ x \pt) + X is given by 

^ = E / fi/pMfi/pM E (oi), (i) 

where fc/ p and fj/p denote the parton densities and A specifies the helicity 
state. At non-vanishing transverse momentum, the leading partonic subpro- 
cesses i+j^ cc[n] + k occur at order a 3 . To obtain the short- distance cross 
section for a given cc state n, we expand, in the rest frame of the heavy quark 
pair, the partonic amplitude in the relative momentum of the heavy quarks 
and decompose the amplitude in spin and color. The amplitude squared for 
can now be written as 

M = Y,M%.,.{0^i), (2) 

n 

where the sum is ordered as an expansion in v 2 . The cartesian tensors 

(oil) = E (o| xM \^ {X)X ) ^ {x)x \ ^<i..x |0> (3) 

x 
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describe the nonperturbative transition of the cc pair into a polarized quarko- 
nium ip and any number of light hadrons X. Up to -^-corrections to the color 
singlet contribution, Lorentz decomposition of these tensors, together with 
spin symmetry, allows us to reduce the matrix elements to four nonpertur- 
bative parameters (0f( 3 5i)), (0?( 3 5i)), (Of^So)), (C^( 3 P )>- They are 
defined in |L| and include a sum over polarizations A. We do not consider 
color singlet operators such as Vf( 3 Si), in the notation of Q, whose contri- 
butions are expected to be irrelevant. 

For unpolarized production, the calculational procedure sketched above 
leads to results identical to those obtained with the standard amplitude pro- 
jection method @. The polarized production cross section, however, requires 
special care as interference contributions from intermediate cc pairs in differ- 



ent angular momentum states JJ Z no longer cancel || [TO], |llf . In this case, 
the relevant projections are unambiguously specified by rotational invariance 
and spin symmetry implemented in the decomposition of (0^ k l ) and can 
be found in ||. For the S = 0, L = intermediate state, each helicity state 
contributes one third of the unpolarized cross section. For the P-wave inter- 
mediate state, one must first sum over all orbital angular momentum states 
L z and then project with the polarization vector e(A) of the quarkonium. 
The short-distance coefficients that enter da\^ [n] can now be written as 

Ayfnl + Byfn] (e(A) • hf + C^n] (e(A) • k 2 f + D y [n] {e(X)-h) (e(X)-k 2 ), (4) 

where k\ and k 2 are the momenta of the initial state partons and ij = 
99-, 91i 9Qy <?<?• Detailed expressions for the coefficients A, . . . ,D will be pre- 
sented elsewhere. Summing over all polarizations A, we find agreement with 
the unpolarized cross sections computed in 0. 

The nonperturbative matrix elements (0%) that enter the color octet con- 
tributions can be adjusted to fit the experimental data as shown in Figs. [I] 
and 0, where we compare the transverse momentum distributions for unpo- 



larized direct J/ip and ip' production with the most recent CDF data [12 
{\fs = 1.8TeV).[] At large p t , the cross section is dominated by gluon frag- 
mentation into color octet 3 Si charm pairs 0, [13|. We have resummed the 



1 For direct J/ip production, indirect contributions from \ decays are removed, but those 
from ip' decays are included in both the data and theoretical prediction. The indirect ip' 
contribution amounts to about 11% p2[ . In accounting for this contribution, we assumed 
that it is independent of pt ■ 
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Figure 1: Fit of color octet contributions to direct J /ip production data from 
CDF (y/s = 1.8 TeV, pseudorapidity cut \r]\ < 0.6). The theoretical curves 
are obtained with CTEQ4L parton distribution functions, the corresponding 
A4 = 235 MeV, factorization scale /1 = (pi + 4m 2 ) 1//2 and m c = 1.5 GeV. The 
fitted color octet matrix elements are given as the central values in Tab [I]. 

leading logarithms (a s lnp 2 /(2m c ) 2 ) Tl in this channel by solving the Altarelli- 
Parisi evolution equation for the corresponding fragmentation function. Evo- 
lution decreases the short-distance cross section and enhances the fitted value 
of (Og ( 3 5i)). The color octet and 3 Pj channels are significant in the re- 
gion p t <> 10 GeV, but fall as da /dpi ~ l/p^ and become negligible for 
the unpolarized cross section at large p t . Because the 1 S , q 8 ' ) and 3 Pj 8 ^ short- 
distance cross sections have a similar p t dependence, the transverse momen- 
tum distribution is sensitive only to a combination 

MtfSjP'PlP) = (OtCSo)) + — 2 (Of( 3 P )). (5) 

As pt increases from 4 to 20 GeV, k varies from 3.9 to 3.0 for m c = 1.5 GeV. 
Since the and 3 Pj 8 ^ channels are significant only at low pt, we choose 
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Figure 2: Same as Fig. [1] for prompt ip' production. 



k = 3.5. [Note that k = 3 was chosen in 0.] The theoretical prediction then 
depends on two phenomenological parameters. 

As in [0, we fit these two parameters to the ^-distribution, but attempt 
to estimate, at least partially, the theoretical uncertainties by varying the 
parton distribution function (PDF) set and renormalization (factorization) 



scale between 1/2 Jpf + Ami an d 2 J pi + 4m;?. The result is shown in Tab. [T] 



for three PDF sets |L4|], where the first error is statistical only and the sec- 
ond reflects the scale variation. In all cases the coupling a s (/i) is evaluated 
according to the A 4 specified by the chosen PDF set. In case of MRS(R2) 
this implies that a s (fi) is evolved with two- loop accuracy. 

While the value of (Og ( 3 Si)) is relatively stable with respect to the un- 
certainty related to the PDF, the fit of M^ 5 ( 1 So 8 ' ) , 3 Po 8 ' ) ) depends sensitively 
on all effects that modify the slope of the pt-distribution, such as the small- x 
behaviour of the gluon distribution or the evolution of the strong coupling. 
A steeper gluon distribution increases the slope of the ^-distribution and 
leads to smaller values of M'^^Sq^^Pq 8 ^). In addition, next-to- leading con- 



5 



CTEQ4L GRV (1994) LO MRS(R2) 



{Oi^CSJ) 1.06 d= 0.14±5;gi 1.12 ±0.1418:1 1.40 ± 0.22±J:f| 

r J/^fl Q (S) 3 p(8)^ /I Qfi 4-11 c;+l-52 Q on -U 1 1/1 +1-46 -| n O 4- H7+2.79 

26 



M^Y^ V^cT) 4.38 ±1.151^ 3.90 ±1.141^ 10.9 ± 2.071?; 

(Of( 3 -Si)) 0.44±0.08lg;g 0.46 ±0.081^ 0.56 ± O.lllg;! 
AffeC 1 ^ , 3 P (8) ) 1-80 ± 0.56±8:io 1.60 ±0.5118;^ 4.36 ±0.961^ 



Table 1: NRQCD matrix elements in 10" 2 GeV 3 . First error statistical, 
second error due to variation of scale. Ratio of ip' to J / ip fixed. 

tributions in a s and also systematic effects inherent to NRQCD, such as the 
inaccurate treatment of energy conservation in the hadronization of the color 
octet quark pairs, modify the shape of the ^-distribution, but are not re- 
flected in the theoretical errors quoted in Tab. [1]. In view of this, the values 
of M3 '^Sq 8 ^, 3 'Pq 8 ^) should be considered with caution. 

For if)' production, an unconstrained fit of (Of ( 3 Si)) and Mf 5 ( l SQ 8 \ 3 PQ 8 ^) 
leads to a value of the second parameter compatible with zero (with large er- 
ror), or slightly negative for GRV and CTEQ4L PDFs. We conclude, that the 
present data does not allow a meaningful fit of this parameter. We therefore 
performed a combined fit of J / ip and ip' data under the (arguably plausible) 
assumption that the ratio M^^S^^Pq 8 ^)/ {O^ ( 3 Si)) is the same for J/ip 
and ip' . The corresponding results for if)' are also shown in Tab. |l]. 

We next consider charmonium polarization in the (s-channel) helicity 
frame ('recoil frame'). In this frame the polarization axis in the ip rest frame 
is defined as the direction of the ip three-momentum in the hadronic cms 
frame. Let P = p p + pp be the sum of the initial hadron four-momenta and 
Q the if) four-momentum. In the ip rest frame = e(A = 0) = (0, —P/\P\). 
The covariant expression reads 

€ l(Q)» = /,„ = — T7 - ~~5 7\ P n > ( 6 ) 



(P-Q) 2 -M 2 s 

where M is the quarkonium mass, taken to be 2m c , and y/s the hadronic 
cms energy. We then obtain the longitudinal polarization fraction £ = 
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Pt/GeV a L b c c L 

4 0.144 0.457 1.81 0.896 

8 0.059 0.147 0.480 0.191 

12 0.031 0.074 0.230 0.085 

16 0.019 0.047 0.144 0.051 

20 0.013 0.033 0.099 0.034 

24 0.010 0.024 0.073 0.025 



Table 2: Coefficients entering the longitudinal polarization fraction. Pa- 
rameter specifications and pseudorapidity cut as for Fig. |l| 

da^ x=o y J2x da^ x \ The polar angular distribution in ip — > / + /~ decay is given 
by " dT 

— oc 1 + a cos 2 9, (7) 

a cos 9 

with a = (1— 3£)/ (1+0 and 9 the angle between the lepton three-momentum 
in the ip rest frame and the polarization axis, the ip direction in the hadronic 
cms frame (lab frame). Note that in this frame scalar products like €l{Q) ■ &i 
[see Eq. (ffl)] depend explicitly on parton momentum fractions Xi,Xj and can 
not be expressed in terms of partonic Mandelstam variables alone. For this 
reason, the partonic helicity amplitudes for the 3 Sx channel computed in [|7| 
can not be used for the present purpose. 

1 In the large-p t limit, the 3 S[ production channel yields transverse polar- 
ization, while the 1 ( Sg an d 3 Pj channels both yield unpolarized quarkonia in 
this limit. However, for the longitudinally polarized cross section k, defined 
in Eq. (^), varies from 5.9 to 3.1 as p t increases from 4 to 20 GeV, so that 
in the important region of low p t a combination of (Of( l So)) and (Of ( 3 -Po)) 
different from M^^Sq^^Pq 8 ^) is required. We represent the longitudinal 
polarization fraction £ as 

f = ^(PtlM + b jM^M + CL(Ot( 3 Po))/mt 
(OtPSi)) + HOt? S )) + c{Ot{*P Q ))lml 

The color singlet contribution is negligible. The coefficients ai,ci,b,c are 
shown in Tab. |2| as function of px- Note that bi = b/3 exactly. Fig. |^ dis- 
plays a, defined in Eq. (^) as function of p t . The shaded band is obtained 
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Figure 3: a as a function of pt- Parameter specifications as for Fig. |l[ 
The band is obtained from Eq. (§) and the 0(af) fragmentation contribu- 
tions from || as described in the text. The solid curve is obtained with 
Mf^S^^P^) = and shows the importance of this parameter for polar- 
ization. The band does not include theoretical uncertainties in the short- 
distance coefficients, Tab. |^. 

as a combination of the uncertainty (statistical only) in the extraction of the 
NRQCD matrix elements [Tab. [I] and the limiting cases that either <C^(%)) 
or (Ot( 3 P )) is set to zero in the combination M^S^ , 3 Pq 8) ). Note that 
a as a function of p t is identical for J/tp and ip' production, because the 
ratios of NRQCD matrix elements for J/ip and ip' have been fixed. Thus, 
in comparing Fig. ^] with direct J j ip polarization data, ip' feed-down must 
be removed in addition to \ feed-down. At low p t ~ 5 GeV, the theoretical 
prediction is compatible with unpolarized quarkonia. As pt increases, the 
angular distribution becomes more and more anisotropic. At pt ~ 15 GeV, 
the fragmentation contributions of order a A s computed in |J become equally 
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important as the mechanisms computed here.Q This second source of lon- 
gitudinal polarization is included in the curves in Fig. || 



In conclusion, we have evaluated the dominant sources of direct charmo- 
nium depolarization in pp collisions at moderate transverse momentum. In 
the lower p t range about 5 GeV, which, due to limitation in statistics, will be 
most relevant in the present Tevatron analysis, an essentially flat angular dis- 
tribution is predicted, unless the (Of ( 1 5'o)) and (Of ( 3 -Po)) matrix elements 
are significantly smaller than their best fit values, obtained from the pt dis- 
tribution of the unpolarized cross section. With increasing p t , the angular 
distribution is predicted to change rapidly due to the onset of substantial 
transverse polarization. Observation of a pattern as in Fig. |3] would provide 
support for the color octet charmonium production picture. 



While this paper was in writing, charmonium polarization at moderate p t 
has also been addressed by Leibovich Unfortunately, contrary to what 



is stated in the paper, he defines the polarization axis as the ip direction in 
the partonic (rather than hadronic) cms frame. His results for the polar an- 
gle distribution can therefore not be used for comparison with Tevatron data. 



We thank M. Mangano and V. Papadimitriou for discussions and for 
providing us with CDF data. 
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